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The near threshold cp meson production in proton-proton and tt p collisions is 
studied with the assumption that the production mechanism is due to the sub-Ncp- 
threshold iV*(1535) resonance. The 7r°, r\ and p°-meson exchanges for proton-proton 
collisions are considered. It is shown that the contribution to the pp — * ppcj) reaction 

O ; 

from the t-channel 7r° meson exchange is dominant. With a significant N*(1535)Ncp 

coupling (<7jv*( 15 35)A^/47r = 0.13), both pp — > ppcj) and ir~p — > ncp data are very well 

£> ' reproduced. The significant coupling of the iV*(1535) resonance to Ncp is compatible 

O ' 

with previous indications of a large ss component in the quark wave function of the 

G\ 

iV*(1535) resonance and may be the real origin of the significant enhancement of the 

in 

<^ ' cp production over the naive OZI-rule predictions. 
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I. INTRODUCTION 

The meson production reaction in nucleon-nucleon collisions near threshold has the potential 
to gain new information on hadron properties and the experimental database on meson 
production in nucleon-nucleon collisions has expanded significantly in recent years. On the 
other hand, the study of the strangeness content of the quark wave functions of baryons and 
baryon resonances, not only in experimental side but also in theoretical side, has been an 
interesting area 2j, which is expected to provide new information on the configuration of 
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baryons and baryon resonances. In the naive quark model, the nucleon and nucleon resonances 
have no strangeness contents, whereas the meson is an ideally mixed pure ss state. From 
the point of view of the naive quark model the pp — > pp<fr reaction involves disconnected quark 
lines and is an Okubo-Zweig-Iizuka (OZI) rule [3j suppressed process. The study of meson 
production in nucleon-nucleon reactions may provide information on the strangeness degrees 
of freedom in the nucleon or nucleon resonances and is of importance both experimentally and 
theoretically. 

Several years ago, the exclusive production cross section for mesonpro duct ion in pp colli- 
sions at Pi a b = 3.67 GeV/c was measured by the DISTO Collaboration [4|, and the preliminary 
result at an excess energy of 18.5 MeV above the threshold was also published by the ANKE 
group [51. With these experimental information about this reaction, several theoretical pa- 

Trull fi 

pers [6|, 13, |8|, |9] by using various models were published to try to explain the experimental data. 
Recently, more data at other energies are available from the ANKE facility 1(3]. Comparing 
the data for the uo meson production from literature, a significant enhancement of <p/uo ratio of 
a factor 8 is found compared to predictions based on the OZI rule. This findings require more 
theoretical work to understand its origin. 

It is well-known that the iV*(1535) resonance couples strongly to the r]N channel. Recently, 
it was found that the iV*(1535) resonance has a significant coupling to KA in the analysis 
of the J/ip — > pAK + decay and the pp — > pAK + reaction near threshold [ill]. The analyses 
[l2, rJ of the recent SAPHIR and CLAS 7p — > K + A data Q, [l5] also show a large coupling 
of the iV*(1535) to KA. In a chiral unitary coupled channel approach it was found that the 
TV* (1535) resonance is dynamically generated as a pole in the second Riemann sheet with its 
mass, width, and branching ratios in fair agreement with experiments and the couplings of 
the iV*(1535) resonance to KT,, r]N and KA are large compared to the ttN channel [16J. The 
analyses of data on the rf photo-production on the proton for photon energies from 1.527 to 
2.227 GeV also suggest the coupling of the r/'N channel to the iV*(1535) resonance 

From the naive quark model, both rj meson and rf meson have a ss component. It seems 
that the iV*(1535) couples strongly to mesons with strangeness or with ss components. These 
phenomena indicate that there may be a significant ss configuration in the quark wave func- 
tion of the iV*(1535) resonance. So, we expect that the iV*(1535) resonance may also have a 
significant coupling to the (pN channel. 

In this paper, we assume that the productions of the meson in proton-proton and n~p col- 
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lisions are predominantly through the excitation and decay of the sub-0iV-threshold A*(1535) 
resonance. By using this picture, we calculate the pp — > pp<fi and ir~p — > ncf) reactions in the 
framework of an effective lagrangian approach. By comparing with the experimental data we 
find that the coupling of the A* (1535) resonance to the cf>N channel needs to be somewhat 
larger than its the coupling to N p channel. The significant coupling of the A* (1535) resonance 
to Ncf) is compatible with previous indications of a large ss component in the quark wave func- 
tion of the iV*(1535) resonance and may be the real origin of the significant enhancement of 
the cf> production over the naive OZI-rule predictions. 

In the next section, we will give the formalism and ingredients in our calculation, then 
numerical results and discussions are given in Sect. 3. A short summary is given in the last 
section. 



II. FORMALISM AND INGREDIENTS 

We study the pp — > pp<p and n~p — ► ncf) reactions near threshold in an effective Lagrangian 
approach. We assume that the near threshold cf> productions in proton-proton and ir~p collisions 
are through the intermediate excitation of the sub-0A-threshold iV*(1535) resonance. The 7r°, 
r) and p°-meson exchanges are considered for proton-proton collisions. The basic Feynman 
diagrams for the pp —>■ ppcf> reaction and the s-channel diagram for the ix~p — > ncf> reaction are 
depicted in Fig. [TJ and Fig. [2j respectively. 

We use the commonly used interaction Lagrangians for irNN, nNN and pNN couplings, 

C-ttNN = -ig-KNNUNlbT 4 ™ N , (1) 
£ V NN = -WttNNUn^Un, (2) 
£ P NN = -gpNNUN(l» + 7^ °pvd U )T ■ p^U N . (3) 

At each vertex a relevant off-shell form factor is used. In our computation, we take the same 
form factors as that used in the well-known Bonn potential model [l8( 

K N (k 2 M ) = (4) 

M K M 

with n=l for tt° and n- meson; n=2 for p°-meson. ku-, wim and Am are the 4-momentum, mass 
and cut-off parameters for the exchanged- meson (M), respectively. The coupling constants and 



the cutoff parameters are taken as [18|, [l9|, [20|: gl NN /4n = 14.4, , g 2 pNN /^ = 0.9, A„. = A 
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FIG. 1: Feynman diagrams for pp —* reaction. The diagram on the left shows the direct process, 
while that on the right shows the exchange one. pi (£=1,2,3,4,5) stands for the 4-momentum of the 
initial and final particle; k and q stand for the 4-momentum of exchange meson and the intermediate 
resonance (iV*(1535)), respectively. 




n: 




q 



FIG. 2: Feynman diagram for tt p — > n<j) reaction. p n , p p , p^, p n and q stand for the 4-momentum 
of 7r~, proton, (j), neutron and intermediate resonance (iV*(1535)), respectively. 

= 1.3 GeV, A p = 1.6 GeV, and k = 6.1. The value of t]NN coupling constant is extremely 
uncertain, with values of g^ NN /4:7i- between and 7 being quoted in the literature, we use 
g^ NN /A7r = 0.4 since many authors say that it is small (see e.g. 21] and [3]). 

To calculate the invariant amplitudes of diagrams in the Fig. [T]and Fig. [2]with the iV*(1535) 
resonance model, we also need to know the interaction Lagrangians involving the iV*(1535) 
resonance. In Ref. 23], a Lorentz covariant orbital-spin (L-S) scheme for N*NM couplings has 
been illustrated in detail. With this scheme, we can easily write the effective N*(1535)Ntt, 
N*(1535)Nr], N*(1535)Np and iV*(1535)iV0 couplings, 



£itNN* — WN*NttUnUn* + h.c, 



(5) 
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£r,NN* = igN*N v UNUN* + h.C, (6) 
*pNN* = WN*NpUNl5{lti ~ q>ll 2 qU )£ fl (Pp)uN* + h.C, (7) 



q- 

£<I>nn* = igN'N^UNlbdp - 9j ^Y^- ) £ ^{P<i>) u N* + h.c. (8) 

Here and u^* are the Rarita-Schwinger spin wave functions for the nucleon and AT* (1535) 
resonance; £^{p p ) and £ M (p</,) are the polarization vectors of the p and 0-meson, respectively. 
It is worth noting that since the spins of the p meson and <fi meson are 1, both S-wave and 
D-wave L-S couplings are possible for the iV*(1535)iVp and iV*(i535)iV0 interactions. It was 
found that the S-wave coupling has a significant contribution to the partial decay width of the 
N* (1535) resonance compared with the D-wave 2J, |25j. In our calculation we only consider 
the S-wave iV*(1535) resonance couplings to Np and neglect the D-wave iV*(1535) resonance 
couplings to N<p for simplicity. The monopole form factors for iV*(1535)-iV-Meson vertexes are 
used, 

F N*N ( ,2 \ _ A *M ~ m M f q\ 

r M \ K M) — \*2 _ 1,2 ' \ y J 
ly M M 

with A; = A; = A; = 1.3 GeV. 

The Af*(1535)A%, A^*(1535)A^?7 and A^*(1535)A r p coupling constants are determined from 
the experimentally observed partial decay widths of the iV*(1535) resonance. With the effective 
interaction Lagrangians described by Eq. 
and r7v*(i535)^Ar^ can be easily calculated 
decay widths, 

3g 2 N * N7T (™ N + EpyT . . 

r^ {15 35)^ - , (10) 

gW^v + K)pT 

J- N*(i535)^Nr] — 7ZTT7 ' y LL ) 



5]) and Eq. ([6]), the partial decay widths rjv*(i535)->7V7r 
24j ] . The coupling constants are related to the partial 



A-nM 



where 



and 



{M 2 N , - (mjy + m^ h ) 2 ){M 2 N , - (m N - m^) 2 ) 

= u ImJ* ' (12) 



For the A^*(1535)A^p coupling constant, we get it from the partial decay width 
Tat* (1535)^ jVp^jVTTTr, and the partial decay width can be evaluated from the total invariant ampli- 
tude M n* (i535)^Np^Nwn of the N* (1535) — > Np — > Nun decay and a three-body phase space 
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integration, 

■Mn*(1535)^N P ^Ntttt = 9 pirnd N* (1535)N pFp (k p )u N (pi, S 1 )'J 5 ('J p — )x 

G'?(k p )(p 2 -p 3 ) 1 ,u N .(q,s N .), (14) 

,-p TT1 [2 d3 Pl m l d3 P2 1 d *P3 1 

rfl^(1535H7Vp^ - 1^^(1535)^^1 ^ ^ ^8 ^ ^ X 

(27r) 4 ,5 4 (g-p 1 -p 2 -p3), (15) 
where G p iu {k p ) is the propagator of the p meson with the form 

GT(*p) = ~»( g P P '/' )• (16) 
p — nip 

Here g and fc p are the 4-momentum of the A/"* (1535) resonance and the intermediate p meson; 
Pi, rrii, and E\ stand for the 4-momentum, mass, and energy of the nucleon; s\ and Sn* the spin 
projection of the nucleon and the iV*(1535) resonance; P2/3 and E 2 /3 stand for the 4-momentum 
and energy of the final two pions, respectively. In our calculation, we use g 2 /Ait = 2.91 as 



the same as that used in Ref. 26]. 



There is no information for the coupling constant of the iV*(1535)iV0 vertex. We determine 
it from the 7t~p — > ncf) reaction. We assume that the near threshold cf> production in it~p 
collisions is through the intermediate excitation of the sub-0iV-threshold iV*(1535) resonance. 
Then, by comparing the theoretical total cross sections of n~p — > ncf) reaction with experimental 
data, we can extract the N*(l535)N<f) coupling constant. 

In Fig [21 we show the s-channel diagram for the n~p — ► ncf) reaction, the intermediate 
excitation is a sub-n^-threshold A/"* (1535) resonance. Following the Feynman rules and with 
the above formula, we can obtain the invariant amplitude A of the n~p — > ncf) reaction, 



A = gN*NngN*N<t>F N *(q 2 )u(p n , S n )^{^^ U )e fl (p ( f > , S<f,)G N * (1535)(q)u(p p , s p ), (17) 



with s n , Sp, Scj), the spin projection of the cf> meson and the nucleon, respectively. The form 



factor for iV*(1535) resonance, FN*{q 2 ), is taken similar as in Refs. 12|, |27] 

A 4 



FN * {q) -A^V{q*-Ml 



with A = 2.0 GeV. G N , (1535) (q) 
written in a Breit-Wigner form 



N* (1535) / 

is the propagator of the iV*(1535) resonance, which can be 

m, 

C (rA _ 7-9 + ^.(1535) , , 

1 ~ M iv*(1535) + ^ M 7V*(1535)riV*(1535)(s) 
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S 1/2 (GeV) 

FIG. 3: Total cross sections vs the cm. energy for ir~p — * ncj) reactions. The experimental data 



are from Ref. 



29]. 



Here T^, (153 5) (s) is the energy dependent total width of the iV*(1535) resonance. According 
to PDG 24j, the dominant decay channels for the iV*(1535) resonance are nN and rjN, so we 
take 

Pttn{s) , ^ P v n{s) 



TV* (1535) I 



r 



p 7TN {M Nt(1535)J 

where p n (ri)N(s) is the following two-body phase space factor 



^ (1535) -^p^(M 2 



(20) 



TV* (1535)' 



(21) 



From the amplitude, we can easily obtain the total cross sections of the 7r~p — > ncf) re- 
action as functions of the excess energies. By adjusting the iV*(1535)iV0 coupling constant, 
we can compare the theoretical results with the experimental data. Theoretical results with 
9% * (1535) N<j)/ / '^ 7r = 0-13 are compared with the experimental data in Fig. [3], we find an excellent 
agreement between our results and the experimental data. Contributions from the u-channel 
N* exchange and p-meson exchange between the pion and the proton are also checked and are 
found to be negligible. 



With experimental mass (1535 MeV), width (150 MeV), branching ratios of the iV*(1535) 24J ] 
and the total cross sections of the n~p — > ncp reaction, we obtain all the coupling constants as 
listed in Table, fl] 
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TABLE I: Relevant iV*(1535) parameters. 



Decay channel Branching ratios Adopted branching 


ratios g 2 : /47r 


Nit 


0.35-0.55 


0.45 


0.033 


Nrj 


0.45-0.60 


0.53 


0.28 


Np -> Ntttt 


0.02 ± 0.01 


0.02 


0.10 


Ncj) 






0.13 



For the pp — > pp<p reaction, the full invariant amplitude in our calculation is composed 
of three parts corresponding to the iV*(1535) resonance production from 7r°, rj and p°-meson 
exchanges, respectively. 

M= ^2 Mi. (22) 

i=n, rj, p 

Each amplitude can be obtained straightforwardly with the effective couplings and following 
the Feynman rules. Here we give explicitly the amplitude M n , as an example, 

Mir = girNNg^N^N^F^ik^F^^ik^F^iq^ip^S^G^X 
U{P4, S4)75(7m - ^^2 U ) G N*(1535)(q)u(Pl, Sl)ti(p 3 , S 3 )j 5 u(p2, S 2 ) 

+ (exchange term with pi <-> p 2 ), (23) 

where is the spin projection of the meson; Si {i = 1, 2, 3, 4) and p, t [i = 1, 2, 3, 4) represent 
the spin projection and 4-momentum of the two initial and two final protons, respectively. 
G n (k n ) is the pion meson propagator, 

= xj^r < 24 > 

The final-state-interaction(FSI) enhancement factor in the 1 5'o di-proton state are taken into 



account by means of the general framework based on the Jost function formalism 



30| with 



\m- l = k T^-, (25) 
k — ia 

where k is the internal momentum of pp subsystem, and the a and (3 are related to the scattering 
parameters via 

<* + P 2 
a = — r = ——5, (26) 

a/j a + p 



with a = -20.5 MeV/c and (3 = 166.7 MeV/c Q (i.e. a = -7.82 fm and r = 2.79 fm) in the 
present study. 

Then the calculations of the differential and total cross sections are straightforward, 

da(pp ^ ppcj)) = -— f > y.\M\ — — —^6 {pi+p 2 -P3-P4-P5), (27) 

Si Sf 

with the flux factor 



(2tt) 5 v / ( Pi -p 2 ) 2 -m4. (2* 



Since the relative phases among different meson exchanges in the amplitude of Eq. (122]) are 
not known, the interference terms are ignored in our concrete calculations. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

With the formalism and ingredients given above, the total cross section versus excess energy 
e for the pp — > pp<ft reaction is calculated by using a Monte Carlo multi-particle phase space 
integration program. It is known that the near-threshold production of the rj meson in pp — > 
ppt] reaction is thought to occur predominantly via the excitation of the iV*(1535) resonance. 
However, the excitation mechanism of the iV*(1535) resonance in proton-proton collisions is 



currently still debated. For example, Batinic et al. 3l| and Nakayama [32j have found that the 



7r and 77-meson exchanges between two protonsplay dominant roles for the excitation of the 



A^*(1535) resonance. However, Gedalin et al. 33] and Faldt and Wilkin [22| have found that 



the p-meson exchange is the dominant excitation mechanism of the iV*(1535) resonance. Here 
the 7T°, rj and p°-meson exchanges for iV*(1535) excitation are all considered. By using the 
formalism and ingredients described in past section we first study the roles of different meson 
exchanges in the pp — > pp<p reaction. Our calculated results are shown in Fig. H] together with 
the experimental data. The double-dotted-dashed, dotted and dashed-dotted curves stand for 
contributions without the pp final-state-interaction (FSI) from ir°, r\ and p°-meson exchanges, 
respectively. A simple summation of them are shown by the dashed line. One can see that the 
contribution from the t-channel 7r° meson exchange is dominant to the pp — > pp<p reaction in 
our model. The p°-meson exchange has a significant contribution to this reaction, while the 
contribution from the 77-meson exchange is negligible. 

From Fig. H] we can see that our theoretical result without the pp FSI agrees well with the 
experimental data at excess energy e = 83.0 MeV. However, at lower excess energies such as: e 
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FIG. 4: Total cross sections vs excess energies 
tion are compared with experimental data 4 



e) for the pp — > ppcf) reaction from present calcula- 
10|. The double dotted-dashed, dotted, dashed-dotted 
and dashed curves stand for contributions from 7r°, r], /9°-meson exchanges and their simple sum, 
respectively. Solid line corresponds to the results with the 1 So pp FSI. 

= 18.5, 34.5 MeV, the calculated total cross sections are lower than the data by a factor of more 
than 4. It is known that the proton-proton FSI plays an important role for the near threshold 
meson production in proton-proton collisions. We also include the effect of the 1 5'o pp FSI by 
using the Jost-function method |30j in our calculation, the results are shown in Fig. H] by the 
solid line which can reproduce the ANKE total cross section data well. 

The momentum, angular distributions of the <fi meson and the p<f> invariant mass spectrum 
for the pp — > ppcj) reaction at excess energy e = 18.5 MeV and 83.0 MeV are also calculated. 
In Fig. [5] we present our calculated results at excess energy e = 18.5 MeV together with 
experimental data from the ANKE group. Differential cross sections as a function of the cm. 
momentum of the outgoing proton are presented in the upper left panel. The upper right panel 
is the angular distribution of the 4> meson in the total proton-proton cm. frame. The dashed 
lines are pure phase space distributions, while, the solid lines are full calculations from our 
model with the So pp FSI enhancement factor. By comparing with the data, we find that 
the pp FSI plays an important role. Our model can explain the experimental data well. In 
the lower part of FIG. [5] the momentum distribution of the meson and the invariant mass 
spectrum of the outgoing proton and the meson are shown. 

In Fig. El we present our calculated differential distributions at excess energy e = 83.0 MeV 
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FIG. 5: Differential cross sections (solid lines) for the pp — > ppcf) reaction at the excess energy e = 
18.5 MeV compared with the ANKE data [5] and phase space distribution (dashed lines). The upper 
left panel is the momentum distribution of the outgoing proton. The upper right panel is the angular 
distribution of the 4> meson in the total cm. frame; The lower left panel is the distribution of the 
cm. momentum of the <j) meson; The lower right panel is the invariant mass spectrum of the outgoing 
proton and the <p meson. 



together with experimental data from the DISTO group. From our calculation we find that 
there is no need to consider the pp FSI at this energy. An excellent agreement between our 
model calculation and the experimental data both in shapes and magnitudes can be achieved 
without taking the pp FSI into account. This is consistent with ANKE findings at e — 75.9 
MeV. The phenomena may suggest that at excess energy about 80 MeV the contribution from 
pp higher partial waves has already overtaken the 1 5'o partial wave as the dominant contribution 
and the FSI becomes not important. 

In our calculation we only include the contribution of the A^*(1535) in the intermediate 
state. In previous calculations @, Q, [s], Is], the 7rp — > <j)N through t-channel p exchange and/or 
sub-threshold nucleon pole contributions are assumed to be dominant. However these contri- 
butions are very sensitive to the choice of off-shell form factors for the t-channel p exchange 
and the gNNcf> couplings and can be reduced by orders of magnitude within the uncertainties 
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FIG. 6: Differential cross sections for the pp — ► ppcj> reaction at the excess energy e = 83 MeV 
compared with the DISTO data Q]. The dashed line reflects pure phase space, while the solid lines, 
which includes the amplitudes but without the pp FSI. 

of these ingredients. Considering the ample evidence for large coupling of the iV*(1535) to the 



strangeness [11 
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34j | and the iV*(1535) resonance is closer than the nucleon pole to 



the <j)N threshold, it is more likely that the iV*(1535) plays dominant role for near threshold 
<p production from up and pp collisions instead of the nucleon pole or the OZI suppressed (frpir 
coupling. Our calculation with the iV*(1535) domination reproduces energy dependence of the 
7i~p — > <pn and pp — > ppcj) cross sections better than previous calculations. The significant cou- 
pling of the iV*(1535) resonance to Ncf) may be the real origin of the significant enhancement 
of the cf> production from np and pp reactions over the naive OZI-rule predictions. This makes 
it difficult to extract the properties of the strangeness in the nucleon from these reactions pro- 
posed by J.Ellis et al 35]. There are also some suggestions [36, for possible existence of 
an Ncf) bound state just below the Nip threshold. However, contribution of such bound state 
with width less than 100 MeV will give a much sharper dropping structure for the 7[~p — > <pn 
cross section at energies near threshold. If such Ncf) bound state does exist, it should have weak 
coupling to ttN and only gives small contribution to the 7r~p — > <pn reaction. 
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IV. CONCLUSIONS 

In this paper, the near threshold <p meson productions in proton-proton and 7r~p collisions are 
studied with an effective Lagrangian approach. We assume that the production mechanism is 
due to the excitation of the sub-iV0-threshold iV*(1535) resonance following 7r°, 77 and p°-meson 
exchanges between two protons. n°NN, r]NN and p°NN coupling constants (except g V NN) 
and form factors are taken from the Bonn potential model. iV*(1535)iV7r , N*(1535)Ni] and 
N*(1535)Np° coupling constants are determined from the partial decay widths of the iV*(1535) 
resonance. The iV*(1535)iV0 coupling constant is deduced from a fit to the experimental total 
cross sections of the 7r~~p — ► ncf) reaction near threshold with the iV*(1535) resonance model. 
We find that the A r *(1535) resonance has a significant coupling to Ncf) (g%*(i535)N<i>/ ^ 7r = 0-13). 

The total reaction cross sections and differential distributions of the near threshold pp — > ppcf> 
reaction are calculated with the iV*(1535) resonance model without adjustable parameter. Our 
theoretical calculation agrees quite well with experiments near threshold. We find that the 
contribution from the t-channel ir° meson exchange is dominant to the pp — > pp<j) reaction. 

The significant coupling of the iV*(1535) resonance to the cf>N channel together with the 
earlier findings of large couplings of the iV*(1535) resonance to the rjN, rj'N and KA channels 



0, y y 



13, 



24J gives a coherent picture that there is a large component of str ang eness in 
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34 



. It also 



the iV*(1535) resonance as expected by various theoretical approaches jll . 
gives a natural explanation for the significant enhancement of the production from np and pp 
reactions over the naive OZI-rule predictions. For a better understanding of the dynamics of 
these reactions, more experimental data at other excess energies with Dalitz plots and angular 
distributions are desired. 
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